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O Abstract 

We calculate the one-loop flavor violating top quark decay t — > egg in the Minimal Supersym- 
metric Standard Model. We discuss the branching ratios obtained with minimal flavor violation, as 
well as with soft-supersymmetry induced general flavor violation. Based on this rate we calculate 
the cross section for the single top quark production via gluon fusion, gg — ► tc, and evaluate its con- 
tribution to the cross section for single top quark production in pp collisions at the Large Hadron 



O 

\Q . Collider. We calculate all contributions coming from the standard model and charged Higgs loops, 



as well as gluino (and neutralino)-up-type squarks, and chargino-down-type squarks loops. Our 
numerical results show that the gluino and the chargino contributions are largest over the whole 
parameter range in the unconstrained Minimal Supersymmetric Standard Model. While in general 
the gluino contributions dominate the cross section, this result depends on the supersymmetric 

•rH , 

flavor violating parameters in the up and down squark sector, the relative mass of the gauginos, 
and whether or not the Grand Unified Theory relationships between gaugino masses are satisfied. 
In the most promising scenarios, the pp — > tc + tc + X cross section at the Large Hadron Collider 
can reach a few hundreds fb. 
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I. INTRODUCTION 



One of the main goals at the CERN Large Hadron Collider (LHC) is to study the produc- 
tion and decay of top quarks. The importance of studying the physics of the top is obvious. 
It is the quark which is closest to the scale of electroweak symmetry breaking and is there- 
fore most sensitive to that scale, and thus to New Physics (NP) beyond the Standard Model 
(SM). One of the important tests of the SM is its predictions for the yield of single tops 
in hadronic collisions. The measurement of single top production cross sections has turned 
out to be a challenging task so far [l| and only upper limits are obtained. For instance, the 
DO experiment, at Tevatron II with integrated luminosity of 230fb _1 , obtained the following 
upper limits on the s (t)-channel processes (as defined below): 6.4 (5.0) pb, at 95% C.L. It 
is expected that increased luminosity and improved methods of analysis will eventually lead 
to the detection of single top events in Tevatron II and subsequently at the LHC. 

Single top production in hadronic machines has been thoroughly discussed within the SM 
where, at lowest order, one has the tree level contributions of s-channel (qq — > tb through W 
exchange), t-channel (ub — > td via W exchange) and gb — > tW with a top quark exchanged. 
In 2I, [J one finds the most recent SM results, which include Next to Leading Order (NLO) 
corrections. These are predicted to be approximately equal to (all the following cross sections 
are in pb), 6.6 (4.1) for a single t (t) production in the s-channel, and 156 (91) for a single 
t (t) production in the t-channel at LHC [2]. The background for single top production in 
the SM was estimated in 4J. 

At the same time, there has been an increased interest in studying forbidden or highly 
suppressed processes as they appear ideal for finding the physics lying beyond the SM. As 
alluded to before, top quark interactions, in particular, might provide a fertile ground to 
searches for NP. It is expected that if NP is associated with the mass generation mechanism, 
it may be more apparent in top quark interactions, rather than in the light fermion sector. 
Along these lines, there have been suggestions that the Flavor Changing Neutral Currents 
(FCNC) single top quark production could be rather sensitive to non-SM couplings such as 
tcV (V = g, 7, Z) and tcH [5J]. The advantage in looking for FCNC processes in top physics 
is that although these exist in the SM, they are minute, leading to tiny, unmeasurable SM 
effects. In general, any measurable FCNC process involving the top will indicate that one 
is witnessing the effects of NP. Note that here we are only interested in processes that are 
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driven by FCNC couplings, which are highly suppressed in the SM by the GIM mechanism. 
Therefore we do not consider NP corrections to SM couplings (like tbW or Zqq) or the 
contributions of new particles (either external or internal), like Z' or W', except those of 
the supersymmetric partners of SM particles. 

FCNC effects in top production contribute to the following single top production processes 
on the partonic level: eg —> t, eg —> tg, cq(q) — > tq(q), qq — > tc and gg — > tc, as well as 
all the above with c — ► u. These subprocesses have been investigated in the presence of 
FCNC effective couplings and in the framework of various NP models js|. 

Of all scenarios of physics beyond the SM, supersymmetry is the most popular. A 
characteristic feature of supersymmetry is that, in addition to the SM FCNC gener- 
ated by the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix, it can provide large 
soft supersymmetry-generated FCNC which would enhance rates and cross sections be- 
yond SM values. The proton collider LHC can produce supersymmetric particles, such as 
squarks and gluinos, with masses up to 3 TeV; as well as potentially lighter ones, such 
as charginos/neutralinos. Flavor-changing interactions appear in supersymmetry in loops 
involving these particles, and thus enhancements in FCNC signals are expected at the LHC. 

Single top quark production generated through FCNC processes has been discussed within 
the effective Lagrangian formalism in a model independent way [6j], as well as in model- 
dependent scenarios [7| . The purpose of this study is to analyze one such class of rare single 
quark FCNC production: the gluon fusion gg — > tc within the framework of low-energy 
supersymmetry. This process was analyzed in where QCD-only loops (loops of gluino 
and squarks), were evaluated in the context of the unconstrained Minimal Supersymmetric 
Standard Model (MSSM). However it is known from analyses of t — > cV that charginos, and 
sometimes neutralinos, can have a large effect on FCNC. Here we first discuss the rare decay 
t — > egg and show it to be larger than t — > eg over most of the parameter space in certain 
cases. Then we perform a complete analysis of gg — > tc in both the constrained MSSM 
(where FCNC decays and cross sections are driven by chargino-down-like squark loops) and 
the unconstrained MSSM (where gluino and neutralino loops contribute as well). We include 
the SM and charged Higgs contributions, contributions from chargino, neutralinos and gluino 
loops, as well as interference effects between SM and non-SM effects, in the context of the 
most general left-left, left-right and right-right intergenerational squark mixings. We also 
address the observability of these channels at LHC. 



3 



Our paper is organized as follows: After a description of the FCNC sources in the uncon- 
strained MSSM (Section II), we present our complete analysis of the branching ratio for the 
top quark t — > egg in MSSM, and compare it to the SM case, where t — > egg was shown to 
be larger than t — > eg Q (Section III). Section IV is devoted to the calculation of the gluon 
fusion cross section gg — > tc, as well as the evaluation of the cross section for pp — > tc + X 
at the LHC through gluon fusion. We include a detailed numerical analysis of the vari- 
ous relative supersymmetric contributions from gluino and chargino loops with or without 
Grand Unified Theory (GUT) mass relations, in addition to a comparison of the constrained 
versus the unconstrained MSSM predictions, as well as observability of these channels. Our 
conclusions and prospects for experimental observations are presented in Section V. 

II. FCNC IN THE UNCONSTRAINED MSSM 

In the unconstrained MSSM there are two sources of flavor violation. The first one arises 
from the different mixing of quarks in the d- and w-sectors in the physical bases, and it is 
described by the CKM matrix (inherited from the SM). In the minimal version of MSSM 
(the constrained MSSM) this is the only source of flavor violation. The second source of 
flavor violation consists of a possible misalignment between the rotations that diagonalize 
the quark and squark sectors, and it is a characteristic of soft supersymmetry breaking. We 
work in the most general version of the model and discuss the constrained version as a limit. 
The superpotential of the MSSM Lagrangian is 

W = iiH l H 2 + Yi 3 E x Ue 3 R + Y^H^d^ + Y^H 2 Q%. (2.1) 

The part of the soft-SUSY-breaking Lagrangian responsible for the non-minimal squark 
family mixing is given by 

+Y^Q l H% + YiA^QJI 1 ^ . (2.2) 

In the above expressions Q is the SU(2) scalar doublet, u, d are the up- and down-quark 
SU(2) singlets (Q,u,d represent scalar quarks), respectively, Y u ^ are the Yukawa couplings 
and i,j are generation indices. The flavor- changing effects come from the non-diagonal 
entries in the bilinear terms M? M? and M? and from the trilinear terms A u and Ad. 
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Here H 1 ' 2 represent two SU (2) Higgs doublets with vacuum expectation values 
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246 GeV, and the angle (3 is defined by tan/5 = V2/V1, the ratio 



of the vacuum expectation values of the two Higgs doublets and fi is the Higgs mixing 
parameter. 

Since we are concerned with top quark physics, we assume that the non-CKM squark 
mixing is significant only for transitions between the squarks of the second and third gener- 
ations. These mixings are expected to be the largest in Grand Unified Models and are also 
experimentally the least constrained. The most stringent bounds on these transitions come 
from b — > 57. In contrast, there exist strong experimental bounds involving the first squark 
generation, based on data from K°-K° and D°-D° mixing 11]. 

It is convenient to specify the squark mass matrices in the super-CKM basis, in which 
the mass matrices of the quark fields are diagonalized by rotating the superfields. Our 
parameterization of the flavor-non-diagonal squark mass matrices for the up- and down- 
type squarks, for the MSSM with real parameters, reads as follows, 
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M 2 Q q + m\ + cos 2(3{T q - Q q s 2 w )m 2 z 



M l { u,c,t } = Ml {uct} + m 2 u ^ + cos2PQ t s 2 w m 2 z , 

M l{d,s, b} = M k{d,sM + m ls,b + cos 2(3Q b s 2 w m 2 z , (2.6) 

•Au,c,t = A n>Cjt — fi cot (3 , 

■Ad,s,b = A d)Syb - /itan/3, 

with m q , T q , Q q the mass, isospin, and electric charge of the quark q, m z the Z-boson mass, 
Sw = sin 9yy and 6\y the electroweak mixing angle. In the above matrices we assumed that 
significant mixing occurs between the second and third generations only. 

We define the dimensionless flavor-changing parameters {^ijd)ab {AB = 
LL, LR, RL, RR) from the flavor off-diagonal elements of the squark mass matrices 
Eqs. (12.41) and (12. 5p in the following way. To simplify the calculation we assume that all 
diagonal entries in (M?) iL , (M?) LiJ , (Ml) RL and {M^) RR and similarly for (M?) AB , are 
set equal to the common value M| USY , and then we normalize the off-diagonal elements to 
M| USY QQ, 
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(2.7) 

The matrix M.\ can further be diagonalized by an additional 6x6 unitary matrix to 
give the up squark mass eigenvalues 

(A<r as = r^rv. (2.8) 

For the down squark mass matrix, we also can define M 2 - as the similar form of Eq. ( 12. 8ft 
with the replacement of (M^) AB (A, B = L,R) by (M 2 b ) AB . Note that while SU(2) L gauge 
invariance implies that (M?) LL = K CKM {M 2 ^) LL K^ CKM , the matrices (M?) LL and (M?) LL 
are correlated. Since the connecting equations are rather complicated and contain several 
unknown parameters, we proceed by including the flavor changing parameters {S^ d )ab as 



independent quantities, while restricting them using previously set bounds [11]. 
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Thus, in the super-CKM basis, there are potentially new sources of flavor-changing neutral 
currents: Char gino- quark- squark couplings, neutralino-quark-squark coupling and gluino- 
quark-squark coupling, which arise from the off-diagonal elements of (M? ^)lr 
and (M? ^)rr- Previous considerations of flavor violating decays [ijj] in the MSSM have 
shown that both up and down squarks contribute significantly. Our analysis shows that this 
is the case here too, and which one is dominant depends on the parameters of the model, 
and in particular on the relative mass hierarchy between the chargino and the gluino. 

In the super-CKM basis, the quark-up squark-gluino (g) interaction is given by 

3 

C U u 9 = ^ 9s T r st [u\ (IVr Pl f «£ - ut (Tuf^ P R f «* + H.c] , (2.9) 

where T r are the SU(3) C generators, Pl,r = (1 = F7s)/2, i = 1,2,3 is the generation index, 
a = 1, . . . , 6 is the scalar quark index, and s, t are color indices. In the gluino interaction, 
the flavor changing effects from soft broken terms M~, M? and A u on the observables are 
introduced through the matrix Ty. 

The relevant Lagrangian terms for the quark-down squark- chargino {xt) interaction are 
given by 

2 3 

A^~x+=E E { Ui ( F « iaS K ckmU Pl xt (r D y a d a - ut [g U ol (KckmU Pr xt (Id)* d a 

a=l i,j=X 

+ u t [U a2 (KckmY?**)^] P R xt (r D )°' +3)a rfa} + H.c. (2.10) 

where the index a refers to chargino mass eigenstates. Y u ^ g are the diagonal up- and down- 
quark Yukawa couplings, and V, U are the usual chargino rotation matrices defined by 
U* M^+V' 1 = diag (m -,+ , m~,+ ). The flavor changing effects arise from both the off-diagonal 
elements in the CKM matrix Kckm and from the soft supersymmetry breaking terms in 

rv 

Finally, the relevant Lagrangian terms for the quark-up squark neutralino (x°) interaction 

are 

4 3 f 4 
= EE{^i3-^tan^P L x°(^ 

n=l i=l * 

- fli-^NVna + iiV^tan^ P R x° n (IV)- u a - u % N n4 P R x° n (T « a j , 

(2.11) 
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where N is the 4x4 rotation matrix which diagonalizes the neutralino mass matrix M%o, 
N* My^N^ 1 = diag(m^o, m^o, m^o, m^o). As in the gluino case, FCNC terms arise only 
from supersymmetric parameters in Tjj. 

Most of the previous a naly ses of FCNC processes in the MSSM concentrated on the mass 
insertion approximation [15]. In this formalism, the (§) terms represent mixing between 
chirality states of different squarks, and it is possible to compute the contributions of the first 
order flavor changing mass insertions perturbatively, if one assumes smallness of the inter- 
generational mixing elements (5's) when compared with the diagonal elements. However, 
when the off-diagonal elements in the squark mass matrix become large, the mass insertion 



approximation is no longer valid 
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131 ] . In the general mass eigenstate formalism, the 



mass matrix in Eq. (12. 8p (and the similar one in the down-sector) is diagonalized and the 
flavor changing parameters enter into our expressions through the matrix Tu^. So, in the 
rare top decays t — > egg, the new flavor changing neutral currents show themselves in both 
gluino-squark-quark and neutralino-squark-quark couplings in the up-type squark loops and 
in the chargino- squark- quark coupling in the down-type squark loops. Therefore here, as in 
our previous work [l0|, we use the general mass eigenstate formalism as described above. 



III. t -> egg VERSUS t eg IN MSSM 

We present here the comparative analysis of the rare two and three body top quark 
decays, t — > egg and t — > eg, closely following the discussion in our earlier paper [9J. There, 
we have shown that, within the SM framework, the branching ratio of t — > egg is about 
two orders of magnitude larger than that of t — > eg in SM, a phenomenon which can be 
dubbed "higher order dominance", and which was revealed e.g., in b and c-physics in the 



past. For the detailed discussion, see [9] and the relevant references therein. Even though 
the branching ratio for t — > egg dominates the one for the two body decay t — *■ eg, it is of the 
order of 10~ 9 and still too small to be detected in collider experiments. Any experimental 
signal for such decay would indicate physics beyond the SM. So, our aim in this section is to 
extend the discussion in to a favorable beyond SM framework in which we would expect 
larger contributions due to extra sources of FCNC - the unconstrained MSSM. Note that 
we include the SM contributions as well in our calculations. 

The one-loop Feynman diagrams contributing to t — > egg in the MSSM are given in a 



8 



set of diagrams Figs, [fl El [31 HJ and [5] in the 't Hooft-Feynman gauge (£ = 1) representing 
gluino, chargino, neutralino, Higgs, and ghost contributions, respectively. 1 We did not 
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FIG. 1: The one-loop gluino contributions to (75 — > tc in the unconstrained MSSM in 
the 't Hooft-Feynman gauge. 

show the SM diagrams here (since they appear in 9J]) but we took them into account in the 
numerical evaluation, for both the decays and the production mode. 

As in 9], we choose to use the 't Hooft-Feynman gauge in which the gluon polarization 
sum is J2x e /l(^' ^) e *X^> ^) = ~9tiv- In order to obey unitarity, this simple choice results in 
the existence of QCD ghost fields whose contributions are shown in Fig. [51 We closely follow 
the method outlined in |9| and references therein for handling the ghost diagrams. 



1 Note that we display the one-loop diagrams for the process gg — > tc. The diagrams for the decay can be 
easily obtained by crossing. 
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FIG. 2: The one-loop chargino contributions to gg — > tc in the unconstrained MSSM in 
the 't Hoof t-Feynman gauge. 
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FIG. 3: The one-loop neutralino contributions to gg — > tc in the unconstrained MSSM 
in the 't Hoof t-Feynman gauge. 

Divergences inherent in the t — > egg calculation arc ultraviolet, infrared, and collincar 
types [9|]. In numerical evaluations, we used the softwares FeynArts, FormCalc, and 
LoopTools [16] to obtain our results. In addition to these, HadCalc [17| is used for de- 
riving the pp process corresponding to the gg fusion discussed in the next section. Using 
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FIG. 4: The one-loop charged Higgs contributions to gg — > tc in the unconstrained 
MSSM in the 't Hoof t-Feynman gauge. 
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FIG. 5: The one-loop QCD ghost contributions to gg — * tc in the unconstrained MSSM 
in the 't Hoof t-Feynman gauge. 

utilities offered by FormCalc, we checked ultraviolet finiteness of our results numerically, 
and introduced phase space cuts to avoid infrared and collinear singularities. 2 



2 These cuts lead to some uncertainties in our results. A more precise approach re quir es full consideration 
of the next-to-leading order corrections to t — * eg, similar to the ones in b decays 18 1. 
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Having mentioned some qualitative features of the decay t — > egg, we do not present here 
most of the analytical intermediate results. We do this since the calculations are lengthy 
and uninspiring. Furthermore, we use well known programs. 3 We have also checked our 
calculations with similar ones, whenever published, as we discuss in the next section. 

We express the matrix element squared \A4\ 2 as a sum over the various contributions. 

n 

These include the SM contribution as given in our previous work |9(. From Figs. HJ El El 
HI and [5J we obtain expressions for the following non-SM terms: the gluino contribution, 
chargino, neutralino, charged Higgs and finally the contribution of the ghosts. 



The results were expressed in terms of Passarino-Veltman functions [19j |. Numerical 
evaluations of these functions have been carried out with LoopTools, which does not require 
reduction of Passarino-Veltman functions to the scalars Aq,Bq,Cq and D . The analytical 



expressions are obtained with the use of FeynCalc 20]. 
The partial width dT for the decay t — > egg is given as 

dV{t^cgg) = -—^{M^d^k^^k^k 



2m t 

spins 



/ t 7 7 \ d^k 2 d^k% d^k^ , . a r (A\ f . . . \ / \ 

d^{kvMMM) = j^p^pao W s^ih-h-h-h), (3.1) 

where ki(k 2 ) is the momentum of the top( charm) quark and k 3 ,k^ the momenta of the 
gluon pair. The volume element can further be expressed as 



/* (^3 ) max ri™2 ) max 
d$ 3 (k i; k 2 ,h,h) = -— / 3 dk° 3 / dk° 2 , (3.2) 



where the limits are 



(A;°) min = Max 



Cuit, — 



2 

(A .0 )m ax = ^+JM (1 _ 2C)) 



f : 



(^) mm = Cm 

(^ 3 °) max = ^(1-2C), (3.3) 
with a = m t — k®. In addition, C is the cutoff parameter, chosen nonzero to avoid infrared 



and collinear singularities 



a. 



For the numerical calculations in the rest of our study we fix 



The complete analytical results can be obtained by contacting one of us (I.T.) 
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C = 0.1, which is large enough to be able to reach the jet energy resolution sensitivity of 
the LHC detector. The results are sensitive to the choice of the C parameter; we find that 
by decreasing C to 0.01, BR{t — > egg) can increase by a factor of 2-4. 

The total decay width of the top quark is taken to be T t = 1.55 GeV. The parameters 
used in our numerical evaluation are given in Table [B 



TABLE I: The parameters used in the numerical calculation. 



a s (m t ) 


a(m t ) 


sin 6w{mt) 


m c (m t ) 


m h (m t ) 


m t (m t ) 


0.106829 


0.007544 


0.22 


0.63 GeV 


2.85 GeV 


174.3 GeV 



The MSSM parameters Msusy ; M 2 , m^o, /i, A, and tan/? are chosen as free for the con- 
strained MSSM and the SUSY-GUT mass relations are assumed. 4 (This is the first scenario 
we consider). Inclusion of the flavor violating parameters 8's among second and third gen- 
eration squarks (the unconstrained MSSM) adds eight more free parameters. Imposing 
SUSY-GUT relations favors a heavy gluino, which decreases the gluino contributions for 
both processes under consideration, t — > cg(g) and gg —>■ tc, and which enhances chargino 
contributions, since the lightest chargino becomes much lighter than gluino. 

As a second scenario we consider the constrained and unconstrained MSSM without 
imposing SUSY-GUT relations. In this case, we run the U(l) gaugino mass parameter Mi 
and the gluino mass Mg separately. 5 Thus the two scenarios we concentrate on are MSSM 
with, and MSSM without, SUSY-GUT relations. 

Given the still large number of parameters in either of these scenarios, the parameter 
space needs to be reduced by making further assumptions. So, for simplicity, we assume 
that the soft SUSY-breaking parameters in the squark sector are set to the common value 
^susy- In addition to this, the trilinear linear terms A u . and A^. are chosen to be real and 
equal to each other and \i is also taken to be real and positive. 

4 The existence of a GUT theory at Planck scale leads to relations among gaugino mass parameters of the 
form 

Mi = (5sw/3cw)M 2 = (5a/3(^ v a s )m,g 

where a and a s are running coupling constants. 

5 We still keep the relation between Mi and Mi , rather than fixing them independently, since this does not 
affect significantly the final results. 
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In the case of flavor violating MSSM, only the mixing between the second and the third 
generations is turned on, and the dimensionless parameters <5's run over as much of the 
interval (0,1) as allowed. 6 The allowed upper limits of <5's are constrained by the requirement 
that m- j > and consistent with the experimental lower bounds (depending on the chosen 
values of M SU sy, A, tan/3, and /i). We assume a lower bound of 96 GeV for all up squark 
masses and 90 GeV for the down squark masses 22j. The Higgs masses are calculated with 



FeynHiggs 



2l| . with the requirement that the lightest neutral Higgs mass is larger than 



221 ]: 96 GeV for the lightest chargino, 



114 GeV. Other experimental bounds included are 
46 GeV the lightest neutralino, and 195 GeV for the gluino. Throughout the paper, only 
m^o and A are fixed globally in the decay and production separately, m^o = 400 GeV and 
A = 620 GeV in the decay process t — ► egg (and t — > eg as well) and m A o = 500 GeV and 
A = 400 GeV, respectively, in the single top production process gg —>■ tc. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
(^lr = (i£ 3 )RL 

FIG. 6: Left panel: Branching ratios of t — * egg and t — > eg decays as functions 
of (<5^)ll with the assumption that GUT relations hold. Right panel: Branching 
ratios as functions of (5^)lr = {5^)rl under the same conditions. The parameters 
are chosen as tan/3 = 10, M SUS y = 300 GeV, M 2 = fi = 200 GeV. 

The rest of the section is devoted to the presentation of our results for the three body 
decay t — > egg and the comparison with the two body channel t — > eg, both within the 
MSSM framework. Since the flavor violating parameters 5's play very important role in 



Even though 6's are allowed to be negative, we run them in the positive region. 
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both decays (both are flavor- violating rare top decay channels), we vary them by keeping 
only a single flavor off-diagonal element non-zero unless otherwise stated. In this section, 
tan (3 = 10 is chosen in all figures except for Fig. [SJ where the dependence of the BR's on 
tan j3 are shown. Furthermore, the common SUSY scale Msusy = 300 GeV; M 2 = 200 GeV, 
and /i = 200 GeV are chosen and fixed globally in this section. Since we are only interested 
in the relative size of the BR(t — > egg) with respect to BR(t — > eg), we consider the scenario 
of MSSM with GUT gaugino mass relations for illustration purposes, and present the case 
without GUT mass relations in one figure at the end of the section, namely Fig. [9j 

Fig. [6] shows the branching ratios of the decays t — > egg and t — > eg as functions of 
(^d)ll on the left panel, and as functions of (5^)lr = {&d )rl on the right panel. Since 
the flavor off-diagonal 5's in the up sector are switched off, these figures show chargino-only 
contributions. As seen from the panels, BR(t — > egg) is almost two orders of magnitude 
larger than BR(t — > eg) in most of the parameter space, and especially for small 8^, up to 
5ff ~ 0.4. As J^'s become larger, BR(t — > eg) increases rapidly and becomes larger than 
BR(t — > egg) for (5 2 d)ll > 0.6 for left-panel and for (S^)lr > 0.8 for the right panel. The 
maximum value reached is around 10~ 7 for non-zero (5^)ll and 10~ 8 for the special case 
{$d)lr = {$d)rl- (Note that t — > eg can get even larger in this part of the phase space). 
These two figures demonstrate explicitly that t — > egg is larger than t — > eg over most of 
the parameter space. We have checked the dependence of BR{t — > egg) and BR{t — > eg) 
on (5^)rr and observed that BR(t — > cp^) remains two orders of magnitude larger than 
BR(t — > cp) for the most part of the interval, while the sensitivity to (5j^)nFi variations is 
not as pronounced as in the (depicted) LL and LR, RL cases. In this case, BR{t — > eg) 
can reach a few times 10~ 9 . 

In Fig. [71 the (5^ 3 )ll and (5^ 3 )r_r dependence of the branching ratios of t — > egg and 
t — > eg decays are shown on the left and right panels, respectively. Since the GUT relations 
are assumed to hold, the gluino mass is rather heavy, about 600 GeV, when M 2 is chosen as 
200 GeV. The two orders of magnitude difference between the BR's for the flavor conserving 
MSSM disappear once we introduce a small flavor violation (~ 0.1) between the second 
and third generations in the up squark sector, which holds for either LL or RR case. The 
branching ratio of t — > eg exceeds that of t — > egg for 5"^ > 0.1. The maximum attainable 
branching ratio for t — > eg is around 10~ 7 , and for t — > egg, 10~ 8 — 10~ 7 which represents 
two orders of magnitude enhancement for t — > egg, and more than 4 orders of magnitude 
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FIG. 7: Left panel: Branching ratios of i -> egg and t — > eg decays as functions 
of with the assumption that GUT relations hold. Right panel: Branching 

ratios as functions of {5^j)rr under the same conditions. The parameters are 
chosen as tan/5 = 10, M SU SY = 300 GeV, M 2 = (jl= 200 GeV. 

enhancement for t — > eg, with respect to the constrained case. The case of (5^ 3 )l_r = (5^ 3 )kl 
is very similar to the case with non-zero (<5^ 3 )ll (left panel) or (<5^ 3 ),r_r (right panel). 

Fig. [H] shows the tan (3 dependence of the decays with zero flavor off-diagonal parameters 
S = for Msusy = 300 GeV, M 2 = fi = 200 GeV. For the decay t -> egg, the SUSY 
contribution comes from the chargino sector in the MSSM (there are no gluino or neutralino 
contributions.) Overall the SM contribution dominates over the MSSM one and the tan/3 
dependence is insignificant, as expected, since the constrained MSSM gives smaller contri- 
butions than the SM to FCNC decays at one-loop level. There is a mild dependence on 
tan/3 for t — > eg decay in the very large tan/? region (> 25). In addition to that, we an- 
alyzed the case with non-zero <5's as well and, for example, for (S^ll = 0.4, we obtain 
BR(t — ► egg) almost two orders of magnitude larger than BR(t — ► eg) in the entire tan/3 
interval considered. 

The last figure of the section, Fig. [HI presents the dependence of the branching ratios 
on the SUSY flavor-violating parameters in the MSSM without SUSY-GUT relations. For 
illustration, we present the (S^)ll dependence of the BR's for the gluino mass nig = 200 
GeV on the left panel, and for m g = 300 GeV on the right panel. The other parameters 
are chosen the same as before, Msusy = 300 GeV, M 2 = ji = 200. As seen from the figure, 
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FIG. 8: The branching ratios of t — ► egg and t ^ eg decays as functions of tan/3 with 
the assumption that GUT relations hold. It is further assumed that all flavor 
off-diagonal parameters S's are zero in both the up and down sectors (constrained 
MSSM) . The other parameters are chosen as M S usy = 300 GeV, M 2 = [i = 200 GeV. 

the relative difference between the decays not only disappears immediately after switching 
{$u)ll on (more precisely, for (S^)ll > 0.01) but also t — > eg exceeds t — > egg with a 
constant factor of 5. This is a gluino dominated case which favors the two-body decay 
t — > eg over the three body decay. The decay t — > eg can get as large as 10~ 5 for m~ g = 200 
GeV and 10~ 6 for m g = 300 GeV. 

From the analysis in this section, it is fair to say that the branching ratio for the three 
body t — > egg decay dominates largely over the one for the two body t — > eg mode for 
the flavor conserving MSSM scenario with SUSY-GUT relations, and remains larger even 
if non-zero flavor off-diagonal parameters in the down squark sector are turned on. Such 
dominance is valid only for relatively small flavor violating parameter in the up squark sector 
{{^W)ll < 0.1). Our results here show that the t — > egg channel gives a larger contribution 
(and may be easier to access) than t — > eg channel over most of the parameter space if the 
flavor violation originated from the down squark sector. 

The predictions of the constrained MSSM (without intergenerational squark mixings) are 
similar to the SM ones. Thus the existence of such SUSY FCNC mixings, directly related 
to the SUSY breaking mechanism, is crucial for the enhancement of the branching ratios. 

Another motivation for considering t — > egg is the issue of single top quark production, 
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FIG. 9: The branching ratios of t — ► egg and i — > eg decays as functions of {5"jj)ll 
without GUT relations for m g = 200 GeV, on the left panel, and nig = 300 GeV, on 
the right panel. The parameters are chosen as tan/? = 10, MgusY = 300 GeV, M2 = 
H = 200 GeV. 



which is one of today's challenging task at colliders. If t — > egg is a promising channel 
with respect to t — > eg 7 , the next question would be what are the consequences of this 
for the single top quark searches at colliders. For this purpose, gg —>■ tc + ic needs to be 
considered. Gluons will become very important and abundant at the LHC, which reaches 
very high energies. Therefore, the rest of the paper is devoted to the consideration of the 
pp — > tc + ic + X cross section at LHC, within the flavor- violating MSSM, by assuming only 
the gluon fusion contribution at partonic level. 



IV. pp -» tc + tc + X AT LHC 

Having discussed the decay mode t —>■ egg and shown that it is a more promising signal 
than t — > eg in the previous section, we consider here the top-charm associated production 
via gluon fusion gg — > tc + ic at the partonic level. Since, at the LHC, TeV or even higher- 
scale energies are going to be probed, gluons inside the proton will become very important. 
This process, as well as other channels involving light quarks, has been considered by Liu et. 



al 



in the unconstrained MSSM driven by SUSY-QCD contributions only. Their results 



7 The observability of t — > egg at LHC will be briefly discussed at the end of the Section [IV 
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show clearly that tc production through gluon fusion is the dominant channel over the ones 
involving light quarks qcf, q, q' = u, c, d, s. For example, 87% of the total hadronic cross 
section a(pp — > tc+ic+X) comes from the partonic channel gg — > tc+ic for (5ff) ll,rr = 0.7 
js]. We agree with their results presented in s| once we make the required modifications to 
the input parameters. 

Here, we present the complete calculation of the hadronic cross section a(pp — *> tc+ic+X) 
at LHC by including all one-loop contributions. In addition to the gluino, the chargino, 
neutralino, and charged Higgs loops as well as the SM part is included. The full set of 
Feynman diagrams contributing to the process at one-loop level through gluino, chargino, 
neutralino, and Higgs loops is given respectively in Figs. (TJ [2J, El and H] in the 't Hooft- 
Feynman gauge. As mentioned in the previous section, we did not display here the SM 



diagrams available in our previous paper 



9J for the t —>■ egg decay case. Note that, as 



mentioned before, working in the 't Hooft-Feynman gauge for this process requires the 
inclusion of QCD ghost diagrams, represented in Fig. [5] 

The partonic level differential cross section for gg —> tc can be expressed as 

1 



da = on o^/o Nou t |A4| 2 ^ 3 , 



\p\L t = {1± ^-^ (4-1) 

where SI3 is the angular volume of the third particle and \f§ is the partonic center of mass 
energy. 8 The matrix squared \Ai\ 2 can be calculated by using the expressions, for t —>■ egg 



by simply using the crossing symmetry (see for example 23|]). Then, the hadronic cross 
section is obtained by convoluting the partonic cross section with the parton distribution 
functions (PDF's), f g / p . So, the total hadronic cross section reads 



f l dC 

Jtio 



/ dt—a(ts,a s {v R )) (4.2) 



where <j(£s, a s {hr)) is the total partonic cross section at the center of mass energy = \f£s 
(y 7 ? is the hadronic center of mass energy) depending on the renormalization scale Here 
£0 defines the production threshold of the process. The parton luminosity is defined as 



dC f 1 dx 

~dj = Jt ~^/p( x ^F)f g /p{^/x,fi F ) (4.3) 



For simplicity we assumed m c zero in our analytical, but not in numerical, estimates. 
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where jjLp is the factorization scale, which is assumed to be equal to the renormalization scale 
fin in our numerical analysis. If one needs to sum over all possible partonic subprocesses 
contributing to the particular final state, there will be a sum over PDF's in Eq. (14. 3p . 

We assume that the top quark in the final state will be reconstructed from events and thus 
it is a physical observable. Of course, to identify the hadronic final state requires making a 
series of cuts on the transverse momentum of the top and charm quarks, the rapidity 77, 
and the jet separation Ai? 34 . The following set is used for the cuts 

p Tc ,p Tt > 15GeV rj c ,rj t <2.B AR u >0A. (4.4) 

Their effect is translated into cuts on the limits of £ in the calculation of the partonic 
cross section. For the transformation of the initial partons to initial hadrons, the program 



HadCalc 



171 ] was used, incorporating the Les Houche Accord Parton Density Function li- 



25]. 



brary (LHAPDF) version 4.2 |24| with the recent data set CTEQ6AB 

For the numerical calculations, we have chosen as input parameter the values m^o = 500 
GeV, A = 400 GeV. In addition the hadronic center of mass energy -^i = 14 GeV is taken 
for the LHC. The factorization and renormalization scales are chosen as the production 
threshold of the process (/ip = fiR = 174.93 GeV). 

We discuss the dependence of the total hadronic cross section of gg — > tc + ic process, 
a(pp — > tc + ic + X), on various MSSM parameters for certain 5 values in scenarios with 
and without GUT relations. Note that for simplicity we assume a common 5 parameter in 
the up and down sector (5fj 3 D ) L L = {^ d )rR = ($u,d)lr = {^u,d)rl and only one (U or 
D) non-zero at a time. At the end of this section we discuss the relative magnitude of the 
contributions coming from the gluino, chargino, and the rest. 

Fig. CEO] shows the Msusy dependence of the total hadronic cross section a(pp — > tc+ic+X) 
for tan/3 = 5, /i = 250 GeV, M 2 = 200 GeV, and m- g = 300 GeV. On the left panel, there 
are two curves, for (S^)ab — 0, and 0.4. The cross section depends very weakly on MgusY 
and, for the constrained MSSM case, the SM contribution is the dominant one. There is 
an enhancement of more than 6 orders of magnitude in the unconstrained MSSM over the 
constrained one, and the cross section can be as large as 15 fb for (5^)ab = 0.4. In the 
down sector, shown on the right panel, the sensitivity of a to Msusy is quite strong, and 
there is an enhancement of about two orders of magnitude in the interval 250 — 1000 GeV. 
There are still 2 — 4 orders of magnitude difference between the constrained MSSM versus 
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FIG. 10: The total hadronic cross section a{pp — > tc + tc + X) via gluon fusion as a 



function of M S usy for tan/3 = 5, p, = 250 GeV, M 2 = 200 GeV, and m~ g 



300 GeV. 



On the left panel, Sff = 0.4 is chosen and compared with the constrained MSSM case 



bjf =0. The same is shown on the right panel for <5|j =0.5. 



10 1 



x io- 1 

+ 

u 

+ io- 2 

io 

B 

b 

IO" 4 

io~ 5 

10"" 



(SfiV = • 
(*u 3 )ab = 0.2 • 
(*u 3 )ab = 0.4 ■ 



200 250 300 350 400 450 500 550 600 
Mg/GeV 





io 2 - 




io 1 - 




10° r 






+ 


10- 1 r 


i u 




+ 


10~ 2 - 










f 


io- 3 - 


ft 




b 






10~ 4 - 




io- 5 - 




io -6 - 





(<*u 3 )ab = 0, Mj = 200 GoV 

(<5u)ab = 0.4, M g = 200 GeV 

(*d )ab = 0.5, Mj = 300 GeV 



15 
tan j3 



20 



25 



30 



FIG. 11: On the left panel, the total hadronic cross section a(pp — > tc + tc + X) via 
gluon fusion as a function of nig for tan/3 = 5, MsuSY = fi = 250 GeV, and M 2 = 200 
GeV at various (5^)ab, A,B = L,R. On the right panel, a(pp — » tc + tc + X) as a 
function of tan/3 at various 5]f D and values. 



the unconstrained MSSM scenarios at (5 2 d)ab = 0.5. The maximum cross section is around 
0.1 fb at around M SU sy ~ 250 GeV. 

In Fig. [HI on the left panel, the total cross section a{pp — > tc + tc + X) is shown as a 
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FIG. 12: On the left panel, the total hadronic cross section a(pp — > tc + ic + X) via 
gluon fusion as a function of fi for tan/5 = 5, MgusY = 250 GeV, M 2 = 200 GeV, and 
rrig = 300 GeV at various (<5y 3 )ABjA, B = L,R. On the right panel, a(pp — ► tc + ic + X) 
as a function of M 2 at various (6^ d )ab with GUT mass relations. 

function of the gluino mass for various (<5{) 3 )ab values. Again the constrained MSSM case 
is dominated by the SM contribution, while for the unconstrained MSSM, o ~ 45 fb for 
m- g = 200 GeV and {5]j)ab — 0.4, which is more than 7 orders of magnitude larger than 
for the case with (S^)ab = 0. On the right panel, the tan/? dependence of the total cross 
section is shown for fj, = 250 GeV, M 2 = 200 GeV, and m g = 200, 300, {6$) A b = 0, 0.4, 
and (S^)ab = 0.5. For very large tan/3 values, the cross section reaches 0.001 fb in the 
constrained MSSM, while in the unconstrained case, for (5^)ab = 0.4 and m g = 200 GeV, 
a cross section of 60 fb is obtained. For (S^)ab = 0.5 and m g = 300 GeV, the cross section 
under these conditions reaches a few fb. 

Fig. [T2] illustrates the /1 (on the left panel) and M 2 (on the right panel) dependences of 
the total hadronic cross section a for representative values of (5^ 3 d )ab- The parameters are 
tan/5 = 5, M SUS y = 250 GeV, M 2 = 200 GeV, and m g = 300 GeV for the left panel, and 
tan (3 = 5, MgusY = A* — 250 GeV, and m g = 300 GeV for the right panel. For non-zero 
(S^j)ab, the cross-section a is not sensitive to \i and remains around 15 fb, but it decreases 
significantly with M 2 in the interval M 2 G [150 — 1000] GeV, if there is a non-zero 5 in 
either the up or down sector. The cross section ranges between 10 fb to 0.1 fb for M 2 = 150 
GeV and 1000 GeV, respectively, for {5^) AB = 0.4, and between 0.1 fb to 0.001 fb for 
(5f) AB = 0.5. 
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TABLE II: Relative contributions to the total cross section a(pp — > tc + ic + X), in 
fb, with and without GUT mass relations. (Sfj)AB = {$d)ab = 0, 0.2, 0.4, A,B = L, R 
is considered. The rest of the parameters are A = 400 GeV, tan/3 = 10, MgusY = 
fj, = 250 GeV and M 2 = 200 GeV. For the case without GUT, m g = 300 GeV is used and 
{$\jd)ab are given in brackets. 



{5f)AB = (S%) AB 





0.2 (No GUT) 


0.4 (No GUT) 


Gluino loop 





1.09 (4.05) 


3.07 (14.13) 


Chargino loop 


1.7710~ 5 


0.0052 


0.034 


The rest 


6.1010~ 6 


0.0025 


0.017 



Before concluding, we comment on the relative contributions of the gluino, chargino, 
and the rest (namely, neutralino, charged Higgs, and SM contributions) to the total cross 
section. In Table ILTj we show the relative contributions to a(pp — > tc + ic + X) from gluino, 
chargino, and the rest, in the MSSM with GUT mass relations (the case without GUT 
mass relations is shown for m g = 300 GeV in brackets if different). For simplicity we set 
(S'u)ab = (^d)ab, A,B = L,R and the values 0, 0.2, and 0.4 are considered. The other 
parameters are A = 400 GeV, tan (3 = 10, M SUS y = /i = 250 GeV and M 2 = 200 GeV. The 
case in which no GUT relations between gaugino masses are imposed corresponds (in our 
analysis) to the case in which the gluino mass is allowed to be smaller. This is the reason why 
only the gluino contributions are enhanced in this scenario. The gluino contributions are also 
dominant for the case in which GUT relations are imposed, and the chargino contributions 
are two orders of magnitude smaller. However, one could envisage a case in which the SUSY 
FCNC parameters )ab are allowed to be large, while the ones in the up sector restricted 
to be small or zero, in which the chargino contribution could be dominant. In the case of 
the constrained MSSM only chargino loops contribute. 

So far, we have considered a common flavor violation in the LL, LR, RL, and RR sectors 
of sfermion matrix. However, one could analyze the relative effects of flavor violation in 
each sector and determine how large their contribution to the cross section could be. If the 
flavor violation comes only from the up sector, then (S^)lr is the most sensitive parameter, 
as shown in Fig. [T3l where we took the parameter values m^o = 500 GeV, /i = Msusy = 400 
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FIG. 13: The total hadronic cross section a(pp — > tc + tc + X) via gluon fusion as 
a function of (S'^)lji (other flavor violating parameters are assumed zero) for 
tan/3 = 30, m A o = 500 GeV, A = 300 GeV, \i = M SUS y = 400 GeV, M 2 = 200 GeV, and 



m-g = 200 GeV . 



GeV, A = 300 GeV, trig = 200 GeV, and tan j3 = 30. Of course this is only one of the best 
possible scenarios. As seen from the figure, the total cross section can be as big as 630 fb 
and the gluino contribution becomes dominant if there exist a large flavor violation in the 
up LR sector between the second and third generations. If the flavor violation comes only 
from the down sector, {5 2 ^)ll gives the largest contributions. In this case the cross section is 
dominated by chargino contribution and can be as large as 0.4 fb with the same parameter 
values. 

Finally, we would like to qualitatively comment on the observability of both the decay 
and the production channels of the top quark considered here. The rare decay t — > egg can 
in general be treated twofold way: one can either treat it inclusively with t — > eg or consider 
it as a separate channel. The former means that t — > egg is taken as QCD-correction to 
t — > eg by assuring that two of three final state jets are collinear so that only two can be 
resolved in the detector. The latter can be a competitive possibility if BR(t — >• egg) is 
significantly larger than that of t — > eg. In here, and in our previous work 9J, we conclude 
that t —>■ egg could be potentially more significant than the two-body decay channel in the 
SM 9| and in some part of the MSSM parameter space. However, this should be taken with 
a dose of caution. Collinearity should be avoided by applying certain cuts. The unphysical 
C-parameter introduced in the phase space here plays an essential role to distinguish t — > eg 
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from t — > egg. Even though in our explorations we considered several values of C, C must 
be taken in the range of jet energy resolution of the upcoming LHC detector. For C < 0.1, 
t — > egg dominates over t — > eg for a larger parameter space, thus availability of better jet 
resolution would give an opportunity to detect t — > egg before t — ^ eg. 

At LHC, predominantly tt pairs are produced. If we consider one of the top quarks 
decay mainly as t — > bW and the other one exotically as it -> egg, then the signal would 
be pp — > tt — > (lv)ggcb (4-jets, a lepton and missing energy), where I = e,/i. 9 For the 
single lepton plus jet topology, it is possible to reconstruct the final state fully, and the 
b-quark can be tagged to obtain a cleaner signal, which introduces extra selection efficiency. 
We assume a(pp — > tt) — 800 pb at LHC and also that the W boson decays leptonically, 
not hadronically. Under this conditions, one can calculate roughly the total expected(raw) 
number of events as 

N = a{pp -»• tt) x BR(i -> bW) x BR(W -> Iv) x L x BR(t -> c##) , 

where L is the integrated luminosity which we take as 100 fb _1 . Therefore we have N = 
800 x 10 3 x 1 x (2/9) x 100 x BR{t -> c^) = (1.77 x 10 7 ) x BR{t -> c^). So, one expects 
around (1.8 x 10 7 ) x BR(t —>■ egg) lepton+4 jets events for an integrated luminosity 100/6 -1 . 
However, counting a total efficiency including trigger efficiency, selection efficiency, as well as 
detector geometrical acceptance, one would approximate a total efficiency around 1% 
Thus, the number reduces to (1.8 x 10 5 ) x BR(t — > egg). So, if the flavor violation comes 
from the down squark sector, then for most part of the parameter space t — > egg would 
dominate over t — ► eg by around two orders of magnitude, but both will remain unobserved 
because lepton+jets events are less than a single event. If the flavor violation comes from 
up-squark sector, then t — > eg dominates and can reach 10~ 5 level for light gluino scenarios, 
which might lead an observable event rate around 1.8. Should the integrated luminosity 
increase at later runs of LHC, one could obtain larger event rates (up to 10 events) if the 
flavor violation is driven by the up squark sector. If the flavor violation comes from the 
down sector of the unconstrained MSSM, the event rate remains below the observable level. 

For the single top production case pp — > tc + ic + X, we already included cuts for the 
transverse momentum and rapidity of the charm and top quarks in the final state, as well 
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9 Considering the W decay hadronically would produce a 6 jet final state, requiring determination of the 
multi-jet trigger threshold. 
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as a lower cut for jet separation. In this case, if assume that the top quark is going to be 
reconstructed in the final state one can predict for example 50,000 events for an integrated 
luminosity 100fb _1 and a(pp — > tc + X) = 0.5 pb. A similar total efficiency consideration 
will going to reduce this further but there exist enough events to find a signal under the best- 
case scenario. Anything beyond the above qualitative discussion about the observability of 
the decay and production channels will be considered in more detail in our future paper 27]. 



V. CONCLUSION 



In this study we analyzed two related issues in top quark physics. In the first part of 
the paper, we concentrated on the comparison of two rare top quark decays, t — > egg versus 
t — ► eg, within the unconstrained MSSM, driven by mixing between the second and third 
generations only. To the best of our knowledge, t —>■ egg decay has been considered only 
in our recent study [9] within the SM framework, where BR(t — > egg) was found to be 
two orders of magnitude larger than BR(t — > eg). However, in the SM, BR(t —>■ egg) 
remains at 10~ 9 level, and thus too small to be detectable. Any experimental signature of 
such channel would require the existence of physics beyond the SM which justifies further 
analyses. Here we studied this decay in the MSSM framework by allowing non-zero flavor 
off-diagonal parameters. Our conclusion of the dominance of the branching ratio of t — > egg 
over t —>■ eg in the SM paper remains mostly valid in the MSSM framework, but now the 
BR's can become as large as 10~ 6 — 10~ 5 . 10 For the cases in which we impose the GUT 
relation between gaugino masses bf and assume a flavor violation in down squark sector, 
the large difference in ratio between the t — > egg and t — *> eg modes disappears only in 
the case of very large intergenerational flavor- violating parameters (close to 1, or to their 
maximally allowed upper values). In that case, t —>■ eg exceeds t —>■ egg. In the case of 
non-zero (S'u)ab, t — > eg dominates t — > egg, except in regions of small flavor violation. 
Once we relax the GUT constraints, there is no longer such a large difference between the 
two and three body decays, as long as a small flavor violation is turned on. Once the 
flavor off-diagonal parameters, {5]j D )ab-, are introduced, the difference in branching ratios 
disappears as the parameters are minute and BR(t — > eg) becomes around 5 times larger 

10 The cutoff parameter C=0.1 is being used. 
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than BR(t — > egg). As expected, if the SUSY-GUT relations hold, both modes cannot 
exceed 10~ 7 level (except for flavor violating parameters near their maximum allowed values 
for t —> eg decay), because the gluino mass is large. Once we relax this condition, both 
t — > egg and t — > eg have branching ratios of the order 1CT 6 — and 10~ 5 , respectively. 

Having shown that the three body rare decay t — > egg is indeed important (comparable 
with, or larger than, the two body decay t — > eg), we carried out a complete calculation of the 
single top-charm associated production at LHC via gluon fusion at partonic level within the 
same scenarios discussed above. This production cross section has been considered before 

n 

including only the SUSY-QCD contributions [8]. We performed a complete analysis by 
including all the electro-weak contributions: the char gino- down-type squark, neutralino-up- 
type squark, charged Higgs, as well as the SM contributions. For simplification, a common 
SUSY FCNC parameter 5 is assumed, {5^ D ) LL = (S^ D ) RR = (Sff D ) LR = {5^ D ) RL (in 
the up and down squark sectors), but most often only one common 5 parameter in either 
sector is allowed to be non-zero each time. We have shown that, in the most promising 
scenarios (if a common SUSY FCNC parameter 5 is assumed), the total hadronic cross 
section a(pp — > tc + ic + X) can become as large as 50 — 60 fb and could reach a few hindered 
fb, especially if we relax the GUT relations between the gaugino masses and assume a flavor 
violation from one sector only (LL, LR, RL, or RR). We have shown that the cross section 
could be as large as 600 — 700 fb if a large flavor violation coming from only (S^)l R is 
allowed. 

The comparative gluino, chargino and other contributions to the process have been esti- 
mated. The gluino contributions dominate over most of the parameter space, when allowing 
flavor violation in both up and down squark mass matrices to be of the same order of mag- 
nitude. However, the chargino contribution is non-negligible and would be dominant in 
either the constrained MSSM, or if the flavor violation was allowed to be much larger in 
the down than in the up squark sector. While the chargino loop is two orders of magnitude 
smaller than the gluino for (5^ 3 )ab = ($d)ab 5 the contribution of the "rest" (neutralino, 
charged Higgs and SM) is around half of the chargino contribution, for all of the SUSY 
FCNC parameters chosen. 

Gluon fusion could be more promising than eg — *> t, qq — > tc, or eg — > gt, which we leave 
for a further study {27]. 

The huge differences in prediction between the constrained and the unconstrained MSSM 
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scenarios make LHC a fertile testing ground for the study of SUSY FCNC processes. Any 
significant rate for the top-charm associated production would be a signal of physics beyond 
the SM, and in particular, of new flavor physics. 
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Note added. After submitting the first version of the present paper to the active 
(hep-ph/0601253,vl) we became aware of an analysis similar to ours |28|]. They analyze the 
QCD SUSY contribution, and take as flavor violation only (S^)ll ^ 0, while we switched 
all flavor violating parameters on between the second and third generations. In our analysis 
a non-zero (5^ 3 )lr parameter is dominant (by one order of magnitude over (5^ 3 )ll), which 
makes a comparison of our results to theirs difficult. 
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